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The  effect  of  soil  moisture  on  chlorine  (CI2 )  deposition  was  examined  in  laboratory  chamber  experiments 
at  high  CI2  exposures  by  measuring  the  concentration  of  chloride  (CD )  in  soil  columns.  Soil  mixtures  with 
varying  amounts  of  clay,  sand,  and  organic  matter  and  with  moisture  contents  up  to  20%  (w/w)  were 
exposed  to  ^3  x  104  ppm  CI2  vapor.  For  low  water  content  soils,  additional  water  increased  the  reaction 
rate  as  evidenced  by  higher  CD  concentration  at  higher  soil  moisture  content.  Results  also  showed  that 
the  presence  of  water  restricted  transport  of  CI2  into  the  soil  columns  and  caused  lower  overall  deposition 
of  Cl2  in  the  top  0.48-cm  layer  of  soil  when  water  filled  r*60%  or  more  of  the  void  space  in  the  column. 
Numerical  solutions  to  partial  differential  equations  of  Fields  law  of  diffusion  and  a  simple  rate  law  for 
Cl2  reaction  corroborated  conclusions  derived  from  the  data.  For  the  soil  mixtures  and  conditions  of 
these  experiments,  moisture  content  that  filled  30-50%  of  the  available  void  space  yielded  the  maximum 
amount  of  Cl2  deposition  in  the  top  0.48  cm  of  soil. 

Published  by  Elsevier  B.V. 


1.  Introduction 

For  toxic  chemical  releases,  consequence  assessment  models 
must  account  for  source  terms,  meteorology,  dispersion  rates,  and 
chemical  degradation  (e.g.,  photolysis  and  deposition)  to  accurately 
predict  a  hazard  zone.  Chlorine  (CI2)  is  a  toxic  industrial  chemical 
that  is  of  concern  to  the  transportation  and  defense  communities 
[1],  but  its  fate  is  still  not  well  characterized.  Deposition  has  been 
mentioned  as  a  potentially  important  factor  in  the  fate  of  high- 
concentration  CI2  plumes  when  model  results  are  compared  with 
observations  from  actual  releases  [2,3],  and  one  modeling  study 
focusing  on  dense  gas  deposition  found  that  deposition  may  be 
important  for  some  conditions  but  additional  experimental  investi¬ 
gations  were  needed  [4],  Experimental  measurements  ofCl2  uptake 
on  aerosol  particles  [5,6],  alfalfa  grass  [7],  and  soil  [8]  indicated 
that  deposition  is  fast.  In  addition,  measurements  of  Cl2  deposition 
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during  outdoor  releases  provided  strong  evidence  that  dry  depo¬ 
sition  is  an  important  factor  in  the  fate  of  a  high-concentration 
CI2  plume  under  calm  conditions  [9],  Thus,  deposition  must  be 
included  in  consequence  assessment  models  to  accurately  predict 
the  hazard  area  resulting  from  a  large-scale  CI2  release. 

Previous  laboratory  [8]  and  field  measurements  [9]  demon¬ 
strated  that  soil  organic  matter  affected  the  CI2  deposition  rate, 
so  a  single  deposition  velocity  would  be  inadequate  to  universally 
predict  Cl2  deposition.  An  empirical  relationship  was  derived  to 
predict  the  deposition  velocity  from  the  fraction  of  soil  organic 
matter  and  the  CI2  exposure,  but  this  was  accomplished  for  soil 
blends  with  constant  moisture  content.  The  field  measurements 
showed  a  positive  correlation  between  the  CI2  deposited  and  the 
moisture  content;  however,  since  there  was  also  a  positive  corre¬ 
lation  between  organic  matter  and  moisture  content,  the  effect  of 
moisture  could  not  be  resolved. 

Water  affects  both  the  chemistry  of  Cl2  deposition  and  gas- 
phase  transport.  CI2  reacts  reversibly  with  water  to  produce 
HOC1  and  OCl~,  which  oxidize  organic  molecules,  providing  an 
irreversible  mechanism  for  CI2  deposition  [10],  HC1  is  also  formed, 
which  can  react  with  carbonate  minerals,  providing  another 
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Table  1 

Soil  constituents  (w/w)  and  density. 


Constituent 

Chemical  formula 

Compost  (%) 

Clay  (%) 

Sand  (%) 

Density3  (g/cm3) 

Organic  content 

N/A 

10.2 

2.3 

0.45 

0.9 

Quartz 

Si02 

84.3 

8.5 

94.6 

2.65 

Clinochlore 

(Mg,Fe)6(Al,Si)4O10(OH)8 

36.9 

2.65 

Kaolinite 

Al2Si205(0H)4 

2.2 

2.65 

Muscovite 

KAl3Si3O10(OH)2 

3.4 

2.5 

2.83 

Kyanite 

Al2Si05 

0.4 

3.59 

Mn-Cordierite 

Mn^AUSisOis) 

7.3 

2.66 

Epidote 

Ca2  Al2  FeSi3  Oq  2  OH 

12.0 

3.44 

Piemontite 

Ca2AlMnSi30i20H 

3.2 

3.49 

Albite 

NaAlSi3Og 

21.5 

2.63 

Anorthite 

CaAl2Si20g 

Sanidine 

KAlSi3Os 

3.7 

2.52 

Phillipsite 

(Na.I<,Ca)(Al,Si)408 

4.9 

2.2 

Calculated  density13  (g/cm3) 

2.22 

2.59 

2.63 

■'  Densities  were  obtained  from  [14], 

b  Calculated  densities  determined  from  mass-weighted  averages  of  constituent  densities. 


mechanism  for  irreversible  deposition.  Thus,  regarding  the  chem¬ 
istry  of  CI2  deposition,  higher  water  content  should  act  to  accelerate 
the  observed  rate  of  reaction.  However,  partitioning  into  an  immo¬ 
bile  phase  (in  this  case,  water)  is  also  known  to  slow  transport 
through  a  porous  bed  of  soil  particles  [11],  The  effective  diffusion 
coefficient  for  vapor  transport  is  reduced  by  a  tortuosity  factor  that 
is  related  to  the  void  space  and  water  content  [12,13], 

Therefore,  the  dependence  of  Cl2  deposition  on  soil  moisture 
content  may  not  be  a  simple  linear  correlation.  Here  we  examine  Cl2 
deposition  into  soils  with  moisture  contents  from  0  to  0.2  (w/w)  and 
with  varying  organic  matter  contents.  We  use  the  chloride  ion  (CD ) 
as  a  tracer  for  Cl2  deposition  to  measure  the  depth  of  penetration  of 
Cl2  into  the  soil  columns  with  the  same  nominal  Cl2  exposure.  We 
expect  to  see  evidence  showing  that  water  increases  the  effective 
reaction  rate  of  Cl2  but  decreases  the  transport  rate.  These  compet¬ 
ing  factors  should  yield  a  soil  moisture  of  maximum  Cl2  deposition, 
below  which  Cl2  deposition  increases  with  increasing  moisture  and 
above  which  Cl2  deposition  decreases  with  increasing  moisture. 

2.  Materials  and  methods 

2.1.  Synthetic  soil  characterization  and  preparation 

Synthetic  soil  blends  were  made  by  mixing  basalt  clay  (Welch 
Tennis  Courts,  Inc.,  Sun  City,  FL)  with  sand  and  compost  (the  latter 
two  purchased  from  a  local  hardware  store,  Home  Depot,  Panama 
City,  FL).  Basalt  clay  and  sand  were  used  as  received.  Compost  was 
sieved  (#8  mesh,  <2.56  mm  particle  size)  to  remove  large  parti¬ 
cles.  Mineral  content  and  mass  fraction  of  organic  matter  of  starting 
materials  were  determined  previously  [8]  and  are  listed  in  Table  1. 

All  soil  constituents  were  dried  in  a  90  °C  oven  overnight  to 
remove  moisture.  Dried  starting  materials  were  mixed  and  then 
water  was  added  back  into  the  mixtures  at  the  indicated  mass  load¬ 
ings  (0-0.2).  Soil  types  are  named  for  the  relative  mass  fractions  of 
the  three  starting  materials  in  this  order:  compost:clay:sand.  Thus, 
soil  type  4:1:1  contains  (by  mass)  4  parts  compost  and  1  part  each  of 
clay  and  sand.  All  samples  were  stored  in  sealed  containers  to  min¬ 
imize  water  loss  prior  to  Cl2  exposure.  Synthetic  soil  blends  were 
packed  into  nominally  5  cm  x  10  cm  (diameter  x  length)  stainless 
steel  columns.  The  columns  were  packed  by  adding  a  small  amount 
of  soil  to  the  column,  compacting  the  soil  layer  with  a  hand  plunger, 
and  repeating  until  the  column  was  full.  The  total  volume  of  a 
packed  soil  column  was  163  cm3,  and  columns  were  weighed  to 
obtain  bulk  soil  density  (Pbuik)- 

Column  void  fractions  (cp)  were  determined  from  measured  bulk 
soil  densities  (Pbuik)  and  estimated  densities  for  soils  with  a  void 


fraction  of  zero  (psoiid )-  Psoiid  was  calculated  from  the  measured 
abundances  of  the  mineral  and  organic  content  of  the  starting 
materials  using  0.9  g/cm3  as  the  density  of  the  organic  content  and 
published  mineral  densities  [14]  (see Table  1 ).  Mass-weighted  aver¬ 
ages  of  the  starting  material  densities  were  used  as  psoiid  for  the  soil 
mixtures,  cp  was  then  calculated  using  Eq.  (1). 

^  =  1  _  Pbuik  (1) 

Psoiid 

Table  2  shows  soil  parameters  (p ,  Pbuik.  Psoiid.  organic  content 
normalized  to  the  soil  column  volume,  and  the  maximum  water 
volume  fraction  (flmax)- 

2.2.  Chlorine  exposure 

Six  soil  samples  were  exposed  simultaneously  in  each  exper¬ 
iment.  The  six  samples  were  either  different  soil  mixtures  with 
the  same  moisture  loading  or  the  same  soil  mixture  with  dif¬ 
ferent  moisture  loadings.  For  experiments  with  the  same  soil, 
one  moisture  loading  was  run  in  duplicate  (i.e„  5  different  water 
contents  were  used).  The  soil  mixtures  are  denoted  in  the  first 
column  of  Table  2  according  to  the  notation  described  in  Sec¬ 
tion  2.1.  Six  packed  soil  columns  were  exposed  to  vapor-phase 
Cl2  (chemical  purity  grade,  Airgas  USA,  LLC,  Atlanta,  GA)  at  ambi¬ 
ent  laboratory  temperature  («*20°C)  in  the  deposition  chamber 
described  previously  [8],  The  assembled  apparatus  was  purged 
with  0.5L/min  zero-grade  compressed  air  (Airgas  USA,  LLC)  for 
one  hour  prior  to  introduction  of  Cl2.  Cl2  and  air  were  introduced 
through  separate  ports  approximately  10  cm  from  the  base  plate. 
The  Cl2  flow  rate  was  initially  set  high  to  quickly  reach  a 
Cl2  concentration  ([Cl2])  of  »*3xl04ppm  in  the  chamber  and 
then  adjusted  to  maintain  the  [Cl2[.  [Cl2]  was  measured  in  real 
time  by  a  UV  absorption  cell  described  previously  [8],  [Cl2]  was 


Table  2 

Soil  Parameters. 


Soil  type 

OCa  (g/cm3) 

0b 

Pbuik  (g/cm3) 

Psoiid  (g/cm3) 

^max 

1:0:0 

0.11 

0.51  ±  0.04 

1.08  ±  0.09 

2.22 

0.27 

4:1:1 

0.088 

0.48  ±  0.04 

1.22  ±  0.10 

2.33 

0.30 

1:1:1 

0.061 

0.43  ±  0.04 

1.41  ±  0.11 

2.47 

0.35 

1:1:4 

0.035 

0.43  ±  0.04 

1.46  ±  0.11 

2.55 

0.36 

1:4:1 

0.053 

0.36  ±  0.05 

1.61  ±  0.13 

2.53 

0.31 

0:1:0 

0.041 

0.32  ±  0.06 

1.71  ±  0.16 

2.59 

0.34c 

a  Organic  content. 

b  Void  fraction  (0)  calculated  using  Eq.  (1). 

c  Soil  mixture  0: 1 :0  with  the  maximum  water  content  (0max)  secreted  water  when 
packed. 
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Table  3 

Description  of  variables. 


Variable 

Description 

Units 

c* 

Cg/Cgo 

dimensionless 

Ce 

Gas-phase  CI2  concentration 

mol/[m3gas] 

Cgo  =  2 

CE(z  =  0) 

mol/[m3gas] 

D  =  10-5 

Diffusion  coefficient 

m2/s 

K=0.01 

Solubility  constant 

(mol/[m3water])/(mol/[m3gas]) 

k 

Rate  constant 

[m^soil  column  ]/mol/s 

k'  =  10-500 

kL2/D 

m3/mol 

L  =  0.1 

Soil  column  length 

m 

R 

Soil  constituent  concentration 

mol/[m3soil  column] 

R* 

RIR  0 

dimensionless 

R0= 50-400 

R(r=0) 

mol/ [m3 soii  column] 

t 

Time 

S 

T=  6 

tD/i.2 

dimensionless 

z 

Depth 

m 

z* 

zIL 

dimensionless 

a 

Air  volume  fraction 

[  m3  gas  ]  /  [  m3  soii  column  ] 

0 

Water  volume  fraction 

[m  water ]/[m  soil  column] 

X 

l+e/c/ff 

dimensionless 

T 

Tortuosity  factor 

dimensionless 

<P 

Void  fraction  (6+ a) 

[  m3  gas  +  m3  water  ]  /  [  tTl3  soil  column  ] 

calculated  using  the  Beer-Lambert  law  with  a  Cl2  absorption 
cross  section,  0330,  of  2.52  x  10~19  cm2/molecule  [15]  and  a  path 
length  of  10  cm.  Cl2  exposure  time  totaled  170  min,  during  which 
[Cl2]  was  above  104ppm  for  120 min.  Exposures  ([Cl2]f)  averaged 
(3.9  ±0.9)  x  106  ppm  x  min. 


2.3.  Sample  analysis 

Sample  analysis  was  described  previously  [8],  Briefly,  the  soil 
columns  were  cored  and  vertically  fractionated  with  a  coring 
assembly  consisting  of  a  stainless  steel  tube  with  slits  cut  to  accept 
a  2.5  cm  radius  spatula  and  a  plunger  to  stop  the  corer  when  flush 
with  the  soil  surface.  The  coring  assembly  had  either  coarse  frac¬ 
tions  (slits  at  0.85, 1 .69, 2.54, 5.08  and  7.62  cm  from  the  soil  surface) 
or  fine  fractions  (slits  at  0.48,  0.95,  1.43,  1.91,  and  2.38  cm  from 
the  soil  surface).  Samples  from  each  fraction  (~1  g)  were  obtained 
and  prepared  for  extraction  (10  mL  deionized  water,  18Mf2/cm) 
and  subsequent  analysis  by  ion  chromatography  [16]  to  measure 
CD  concentrations.  Ion  chromatography  was  performed  in  dupli¬ 
cate  on  an  anion  exchange  column  (4  x  250  mm  IonPac®  AS9-HC, 
Dionex®,  Bannockburn,  IL)  with  a  mobile  phase  of  l.OmL/min 
9.0  mM  Na2C03  (pump  model  GP40,  Dionex),  and  CD  was  detected 
by  a  change  in  conductivity  (detector  model  ED50,  Dionex).  A  small 
flow  of  25  mM  H2S04  was  used  to  regenerate  the  suppressor  (model 
ASRS®-ULTRA  II  4  mm,  Dionex).  Standard  solutions  of  NaCl  (99%, 
Sigma-Aidrich)  (0.0178-1.7836  mM)  were  analyzed  to  generate  a 
calibration  curve,  and  injections  of  a  standard  solution  were  made 
after  every  10  samples. 


2.4.  Mathematical  modeling 


Deposition  of  Cl2  into  the  soil  columns  was  modeled  using  Fields 
law  of  diffusion  (Eq.  (2))  and  a  simple  rate  law  for  Cl2  reaction 
with  reactive  soil  constituents  (Eq.  (3)),  which  are  primarily  organic 
matter.  All  variables  are  defined  in  Table  3.  The  quantity  (1  +9l<la ) 
in  Eq.  (2)  is  the  retardation  factor  [11]  (7.)  and  accounts  for  the 
slower  transport  of  a  chemical  species  through  a  porous  material 
due  to  dissolution  into  the  aqueous  phase.  The  tortuosity  factor  (t), 
defined  in  Eq.  (4),  is  an  empirical  correction  to  gas-phase  diffusion 
in  a  porous  material  [12,13], 


8^ 

dt 


tD 


^-K9KCsR 
3 z2  ce  K 


(2) 


'  =  -kOKCoR 
dt  B 


(3) 


(4) 


It  is  helpful  to  use  non-dimensionalized  variables  when  solving 
transport  equations,  so  the  following  variables  were  substituted 
into  Eqs.  (2)  and  (3):  C*,  R*,  T,  and  z*.  In  addition,  k'  was  used,  and 
substituting  these  variables  into  Eqs.  (2)  and  (3)  yield  Eqs.  (5)  and 

(6). 


32r*  9 

r-2—  -k'-KC*R0R* 
3z*2  a 


(5) 


3  R* 
3  T 


=  -9k'I<Cg0C*R* 


(6) 


Average  Depth,  cm 


Fig.  1.  Deposition  of  CI2  (measured  as  [Cl-])  as  a  function  of  depth  using  coarse 
vertical  fractions.  Soil  types,  0,  and  water  mass  loadings  are  indicated.  [Cl-]  is  nor¬ 
malized  to  the  soil  column  volume,  and  error  bars  are  the  standard  deviations  of 
replicate  measurements.  Insets  to  the  top  two  panels  show  a  close-up  of  data  at 
3.8  cm  deep  to  illustrate  the  higher  [Cl-]  detected  in  these  dry  soils. 
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Eqs.  (5)  and  (6)  were  solved  numerically  using  COMSOL 
Multiphysics®,  a  finite  element  solver.  Boundary  conditions  were  a 
constant  CI2  concentration  at  the  soil  surface  (Cgo)  and  zero  CI2  flux 
at  the  base  of  the  soil  column  (z*  =  1 ).  Values  or  ranges  for  variables 
are  disclosed  in  Tables  2  and  3 ;  grid  spacing  ( A z*)  was  0.002.  Values 
for  D  [  1 7]  and  I<  [  1 8 ]  were  obtained  from  the  literature.  U  Cg0,  and  T 
were  set  to  values  used  in  the  experiments.  R0  was  estimated  from 
the  experimental  results,  k'  was  estimated  as  described  in  Section 
3.2,  and  ranges  for  a  and  6  were  obtained  as  described  in  Section 
2.1. 

3.  Results  and  discussion 

3.1.  Chamber  measurements 

Fig.  1  shows  the  effects  of  moisture  content  on  three  soil  types 
analyzed  with  the  coarse  vertical  fractions.  All  soils  exhibit  Cl-  for 
all  moisture  contents  examined.  Cl-  present  in  the  dry  samples 


suggests  that  there  is  some  residual  water  in  the  soil  columns, 
which  is  not  surprising  since  it  is  unlikely  that  all  the  water  is 
removed  from  the  soil  mixtures  at  90  °C.  Significantly  more  CI_  is 
observed  in  all  samples  when  water  was  added.  For  the  1:0:0  soil, 
just  2%  (w/w)  water  loading  yielded  the  same  [Cl-  ]  as  that  with  20% 
(w/w)  water.  This  increased  deposition  of  Cl2  with  added  water  is 
consistent  with  the  chemical  mechanism  of  dissolution  followed 
by  reaction. 

While  the  1 :0:0  soil  exhibits  only  minor  changes  in  the  deposi¬ 
tion  of  Cl2  for  water  contents  from  2  to  20%,  deposition  in  the  other 
two  soil  mixtures  is  strongly  dependent  on  the  water  loading.  Both 
the  1:1:1  and  0:1:0  soils  show  increasing  CI2  uptake  to  8%  water, 
but  decreasing  deposition  at  higher  water  loadings.  This  reduced 
deposition  is  evident  in  decreased  [Cl  ]  in  the  top  fraction  and  those 
deeper.  For  example,  the  1:1:1  soil  with  1 2%  water  content  shows 
a  slightly  reduced  [Cl  ]  in  the  top  fraction  relative  to  that  with  8% 
water,  but  at  1 .25  cm  deep,  there  is  a  significant  reduction  in  the  1 2% 
water  sample.  Further,  the  insets  to  the  1:0:0  and  1:1:1  soil  mixture 


Fig.  2.  Deposition  of  CI2  as  a  function  of  depth  with  finer  vertical  fractions.  Soil  types,  0,  and  0/0  are  indicated.  Water  mass  loadings  are  0, 4,  8, 12,  and  16%  for  ■,  □,  •,  O. 
and  ▲,  respectively.  [Cl-]  is  normalized  to  the  soil  column  volume,  and  error  bars  are  the  standard  deviations  of  replicate  measurements. 
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graphs  in  Fig.  1  show  that  at  3.8  cm  deep  the  dry  sample  contains 
significantly  more  Cl-  than  those  with  added  water.  These  observa¬ 
tions  confirm  that  moisture  does  limit  CI2  transport  through  soils  in 
conditions  similar  to  those  encountered  at  a  hazardous  Cl2  release. 

In  Fig.  1 ,  very  little  Cl-  was  detected  at  depths  larger  than  k.2  cm, 
so  to  provide  a  more  detailed  picture  of  Cl2  deposition  in  soil, 
samples  were  vertically  fractionated  with  finer  resolution  and  addi¬ 
tional  mixtures  were  examined.  Fig.  2  shows  the  results  for  the 
indicated  soil  mixtures  to  a  depth  of  2. 2  cm.  Soil  type  4:1:1,  similar 
to  1:0:0,  shows  no  reduction  in  [Cl  ]  at  the  highest  water  content 
investigated  (16%,  0/0  =  0.48);  however,  at  the  highest  water  con¬ 
tent,  Cl2  does  show  lower  deposition  at  the  deeper  fractions  relative 
to  those  samples  with  less  water. 

While  soil  4:1:1  exhibits  similar  Cl2  deposition  for  moisture 
contents  higher  than  the  dry  case,  the  other  soil  mixtures  in  Fig.  2 
show  significant  effects  of  water  on  Cl2  deposition.  Soils  1:1:1, 
1:1:4,  and  1:4:1  show  a  maximum  deposition  in  the  top  fraction 
at  8%  (w/w)  water.  At  16%  water,  these  soil  mixtures  exhibit  very 
little  penetration.  From  Fig.  1,  Cl2  deposits  to  a  maximum  depth 
of  0.85  cm  for  soil  1:1:1,  and  results  from  the  more  finely  resolved 
fractions  in  Fig.  2  indicate  that  very  little  Cl2  deposits  deeper  than 
0.48  cm.  Clearly,  the  higher  water  content  affects  Cl2  transport.  The 
legends  in  Fig.  2  show  0/0,  which  is  a  measure  of  the  fractional 
pore  space  filled  with  water  and,  therefore,  unavailable  for  gas- 
phase  transport.  Soils  1:1:1  and  1:1:4  show  very  little  deposition 
in  fractions  deeper  than  0.48  for  0/0  =  0.62  and  0.64,  whereas  soil 
1:4:1  shows  significant  deposition  in  deeper  fractions  for  0/0  =  0.61 . 
However,  at  0/0  =  0.85,  Cl2  deposition  in  the  deeper  fractions  of  soil 
1:4:1  is  zero. 

It  is  clear  that  the  void  space  in  the  soil  column  plays  an  impor¬ 
tant  role  in  the  overall  deposition  of  Cl2,  and  the  amount  of  water 
affects  a.  by  filling  the  void  space.  However,  the  water  may  not 
only  fill  pore  space,  it  could  also  affect  the  packing  of  the  soil 
columns  since  samples  were  prepared  with  the  prescribed  amount 
of  water  already  mixed  in.  To  test  whether  better  packing  affects  the 
observed  depth  penetration  of  Cl2,  six  identical  soil  samples  were 
prepared  and  three  of  them  were  dried  in  an  oven  at  60  °C  for  two 
days.  All  six  samples  were  exposed  to  Cl2  at  the  same  time,  so  the 
only  difference  was  the  presence  or  absence  of  water.  This  experi¬ 
ment  was  repeated  for  soil  1:1:1  with  8%  and  16%  water  loadings, 
and  the  results  are  shown  in  Fig.  3. 

For  the  wet  soils  in  Fig.  3,  the  depth  profile  of  [Cl-]  was  simi¬ 
lar  to  that  observed  in  Fig.  2.  Soil  with  16%  water  loading  showed 
very  little  Cl-  in  the  second  fraction  down,  and  soil  with  8%  water 
showed  a  nearly  linear  dependence  of  [Cl-  ]  on  depth.  The  dried  soil 
columns  exhibited  only  a  small  [Cl  ]  in  the  top  fraction,  consistent 
with  the  dry  samples  in  Fig.  2.  In  addition,  [Cl-]  shows  essentially 
no  dependence  on  depth  for  the  dried  samples  in  Fig.  3.  If  water 
significantly  affected  the  packing  of  the  soil  columns,  the  dried 
samples  would  show  a  similar  depth  profile  of  [Cl-]  as  the  wet 
ones.  Instead,  Cl2  transport  is  unhindered  to  a  depth  of  2.2  cm. 
Thus,  the  deeper  relative  penetration  of  Cl2  in  the  dry  samples 
of  Figs.  1  and  2  is  not  due  to  an  artifact  of  the  column  packing 
procedure.  Further,  the  flat  depth  profile  of  the  18%  water  soil 
that  was  dried  in  Fig.  3  demonstrates  that  its  packing  does  not 
hinder  Cl2  transport,  so  the  sharp  drop-off  of  [Cl-]  at  a  depth 
of  0.7  cm  for  the  highest  water  loading  of  soils  1:1:1,  1:1:4,  and 
1:4:1  in  Fig.  2  is  due  primarily  to  water  filling  available  void 
space. 

3.2.  Numerical  modeling 

There  are  several  variables  that  affect  Cl2  deposition,  Cgo,  t, 
R0,  0,  and  0.  In  these  experiments,  Cgo  and  t  were  kept  constant 
but  the  latter  three  parameters  were  varied.  To  corroborate  the 
conclusion  that  water  is  decreasing  the  rate  of  transport  of  Cl2 


Fig.  3.  Deposition  of  Cl2  on  soil  1:1:1  with  (A)  16%  (w/w)  or  (B)  8%  (w/w)  water 
content.  All  soil  columns  were  packed  while  wet,  then  three  columns  in  each  exper¬ 
iment  were  dried  at  60  C  for  two  days  (■)  and  the  other  three  were  covered  to 
retain  moisture  until  the  experiment  (•).  Columns  in  (A)  and  (B)  were  exposed  to 
3.4  x  106  and  3.8  x  106  ppm  x  min  Cl2,  respectively. 


through  the  soil  column,  Eqs.  (5)  and  (6)  were  solved  numerically 
using  COMSOL.  These  equations  use  the  non-dimensional  vari¬ 
ables  T  and  z*,  C*.  and  R*,  defined  in  Table  3.  To  generate  model 
results  that  reasonably  represent  the  conditions  of  these  experi¬ 
ments,  Rq  and  k'  must  be  known.  Extrapolating  the  linear  depth 
dependence  of  [Cl-]  for  soil  1:0:0  (which  contained  the  highest 
amount  of  organic  matter)  in  Fig.  1  to  z  =  0  yields  approximately 
0.28 mmol/cm3,  so  an  upper  limit  of  0.4 mmol/cm3  (400 mol/m3) 
was  used  for  R0.  Soils  0:1:0  and  1:1:4  (both  with  low  amounts  of 
organic  matter)  exhibited  maximum  [Cl-]  at  about  0.1  mmol/cm3, 
so  a  lower  limit  of  0.05  mmol/cm3  (50 mol/m3)  was  used  for  R0. 
These  lower  and  upper  limits  of  R0  adequately  cover  the  experi¬ 
mental  parameter  space.  The  effective  rate  constant  (k')  was  not 
known,  thus  k '  was  adjusted  to  find  agreement  with  the  experi¬ 
mental  results. 

Fig.  4  shows  the  results  of  the  numerical  solution  to  these 
partial  differential  equations,  in  which  the  effects  of  k',  R0,  0, 
and  0  are  shown.  Note  that  the  y  axis  is  the  amount  of  product 
formed  normalized  to  the  product  formed  at  T=  00,  and  the  x  axis 
is  the  normalized  depth.  While  the  focus  of  this  work  is  on  the 
effects  of  water  on  Cl2  deposition,  it  is  important  to  examine  the 
effects  of  other  variables  in  Eqs.  (5)  and  (6)  to  garner  a  compre¬ 
hensive  understanding  of  the  factors  that  affect  deposition  in  soil 
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Fig.  4.  Numerical  modeling  results  for  the  variable  quantities  listed  in  each  quadrant.  (A),  (B),  (C),  and  (D)  show  the  results  for  independently  varying  k',  R0,  <j>,  and  0/0, 
respectively.  Inset  to  (D)  shows  the  integrated  numerical  results  to  illustrate  the  effect  of  sampling  with  Az*  =  0.048.  Other  parameters  were  Cg0  =  2  mol/m3,  T=6,  and  ft  =  0.01 . 


matrices.  In  panel  A,  faster  k'  increases  the  amount  of  product 
near  the  soil  surface  ( z *  <  0.05),  but  decreases  the  amount  of  prod¬ 
uct  at  larger  depths.  The  effective  aqueous-phase  rate  constant 
(Id)  was  unknown,  but  a  value  of  100m3/mol  approximates  the 
depth  dependence  of  the  product  concentration  observed  in  the 
data  of  Fig.  2  with  similar  9/0.  Thus,  k'  =  100  m3/mol  was  used  for 
the  numerical  results  in  the  other  panels  of  Fig.  4.  Panel  B  shows 
that  higher  R0  does  not  affect  the  relative  amount  of  product  at  the 
surface  (z=0),  but  it  does  reduce  the  relative  penetration  of  CI2. 
Panel  C  shows  the  effect  of  void  space  ((/>),  and  a  larger  amount 
of  void  space  yields  more  product  and  deeper  penetration.  Note 
that  in  Panel  C  that  6/0  is  kept  constant  for  the  three  values  of 
<p. 

Finally,  Fig.  4D  shows  the  effects  of  water  on  the  deposi¬ 
tion  of  Cl2.  Increasing  the  water  fraction  decreases  Cl2  depth 
of  penetration  and  increases  the  product  concentration  at  the 
surface  (z  =  0).  As  a  more  direct  comparison  with  the  experimen¬ 
tal  results,  the  inset  to  Fig.  4D  shows  the  numerical  modeling 
results  integrated  over  the  depth  fractions  used  in  the  experi¬ 
mental  analysis.  The  same  trend  is  observed  in  the  numerical 
results  as  those  in  Figs.  1  and  2.  Fligh  water  content  attenuates 
the  transport  of  Cl2  into  deeper  soil  fractions,  and  the  top  frac¬ 
tion  exhibits  significantly  less  product  at  high  water  contents. 
These  results  corroborate  the  conclusion  from  the  experimental 
results. 


As  further  comparison  between  numerical  modeling  and  exper¬ 
imental  results,  Fig.  5  shows  results  from  the  top  0.48-cm  fraction 
for  each.  Experimental  results  (Fig.  5B)  show  peaks  in  the  prod¬ 
uct  yield  at  610  =  03-0.4,  and  the  modeling  results  (Fig.  5A) 
show  peaks  in  the  product  yield  at  progressively  higher  9/0  for 
higher  Rq.  The  9/0  range  for  the  modeling  results  (0.4-0.55)  is 
near  those  of  the  experiments.  However,  the  experimental  results 
show  no  apparent  trend  of  the  910  value  for  the  product  peak 
with  the  organic  concentration  (note  organic  content  in  Table  2). 
This  may  indicate  that  other  processes  or  factors  are  not  ade¬ 
quately  accounted  for  in  Eqs.  (5)  and  (6).  Nevertheless,  the  similar 
trends  observed  in  the  modeling  and  experimental  results  support 
the  conclusion  that  water  restricts  transport  of  Cl2  into  the  soil 
column. 

As  mentioned  in  Section  1 ,  water’s  restriction  of  vapor  trans¬ 
port  through  a  porous  soil  matrix  is  established  in  the  literature. 
However,  much  of  the  environmental  literature  is  concerned  with 
longer  time  scales  and  lower  concentrations  than  that  relevant  to 
the  transport  and  evolution  of  hazardous  chemical  plumes.  The 
experiments  and  numerical  solutions  here  show  that  in  the  case 
of  Cl2,  the  same  mechanism  is  operable  for  deposition  of  Cl2  into  a 
soil  matrix  for  conditions  that  are  relevant  for  hazard  prediction 
and  analysis  (high  [Cl2]  for  a  couple  of  hours).  However,  water 
in  these  experiments  was  stationary  (i.e.,  there  was  no  advective 
transport  of  water).  For  conditions  of  water  flow  (e.g.,  rainwater 
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Fig.  5.  Comparison  of  top  fraction  results  ( Az= 0.48  cm)  from  (A)  numerical  mod¬ 
eling  and  (B)  experimental  results.  Numerical  solutions  were  obtained  using  T=  6, 
0  =  0.5,  and  the  indicated  R0s  (mol/m3).  Experimental  results  are  those  from  Fig.  2 
plotted  against  B/tp. 

seeping  through  the  soil),  the  added  transport  term  would  change 
the  results. 

4.  Conclusions 

CI2  deposition  was  measured  on  synthetic  soil  blends  with  the 
aim  of  elucidating  the  role  of  soil  moisture  content  on  the  rate  of  Cl2 
removal.  Experimental  results  demonstrated  that  water  increased 
the  rate  of  reaction  but  it  decreased  the  vapor  transport  through  the 
soil  column.  The  maximum  Cl2  deposited  corresponded  to  condi¬ 
tions  for  which  the  water  volume  fraction  was  approximately  half  of 
the  total  void  space  in  the  soil  matrix.  Numerical  solutions  to  trans¬ 
port  and  reaction  equations  that  describe  the  problem  confirmed 


the  conclusions.  On  a  final  note,  we  remind  the  reader  that  these 
results  represent  conditions  for  which  condensed  phase  water  is 
stationary.  For  conditions  where  water  is  moving,  the  results  will 
be  different. 
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